Coherent bimolecular reactions with quantum-degenerate matter waves 
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We demonstrate theoretically that the abstraction reaction A + B2 AB + B can be driven coher- 
ently and efficiently with quantum-degenerate bosonic or fermionic matter waves. We show that the 
initial stages of the reaction are dominated by quantum fluctuations, resulting in the appearance 
of macroscopic non-classical correlations in the final atomic and molecular fields. The dynamics 
associated with the creation of bosonic and of fermionic dimer-atom pairs are also compared. This 
study opens up a promising new regime of quantum degenerate matter-wave chemistry. 
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I. INTRODUCTION 

The making and probing of ultracold molecular gases 
have attracted much attention in recent years [ij , as their 
realization opens up exciting applications in fields rang- 
ing from condensed matter physics to quantum infor- 
mation science. The application of magnetic Feshbach 
resonances (FR^ 12, S M @ and of optical photoas- 
sociation (PA) [7|, |8|, |9|, [lOj, oftentimes in combination, 
within an atomic Bose-Einstein condensate can result in 
the creation of diatomic [ll| as well as of more complex 
molecules, as evidenced by recent experimental observa- 
tions of transient Efimov trimcr states Css [l^ and of 
the molecular tetramer CS4 (indirectly via resonances in 
inelastic processes) [13, [13 . We also mention the impor- 
tant work of Ling et al. |l5l.[l6l.[l7t. who proposed to use 
a chirped coupling field to compensate the effects of non- 
linear collisions within the stimulated Raman adiabatic 
passage (STIRAP) technique ^ [U, [H, [H, [M? ? ] 
and thus to efficiently generate large amounts of deeply 
bound ultracold molecules. This is an extension of previ- 
ous developments in associative STIRAP (see, e.g. early 
work by Mackie and coworkers [l^l), that is particularly 
relevant in the context of the present paper. The atom- 
molecule dark states involved in such a process was first 
realized experimentally via coherent two-color PA psj in 
an early example of what is now called superchemistry 

Over the past ten years, ultracold matter- wave su- 
perchemistry has focused on coherent association or dis- 
sociation reactions [13, [1^ between atoms and diatomic 
molecules. For example, Moore and Vardi [131 stud- 
ied the possibility of an almost complete Bose-enhanccd 
channel selectivity in the coherent photodissociation of 
bosonic triatomic ABC molecules, resulting from the in- 
terplay between Bose enhancement and competition be- 
tween modes for a finite number of initial molecules. 
Other recent examples include work on the dependence 
of dissociation on the size or shape of the reaction ves- 
sel (confinement effect) [1^ and the assembly of Fermi- 
degenerate dimers via cooperative association [soj . 



The present paper extends the toolbox of superchem- 
istry to the coherent abstraction reaction (or bimolecu- 
lar reactive scattering) A + B2 ^ AB + B, where A, B, 
B2 and AB denote either bosonic or fermionic atoms or 
dimers. This reaction is an important benchmark sys- 
tem in chemical physics. Its dynamics has attracted 
much interest in recent studies of reactive resonance or 
low-energy non-Born-Oppcnheimer reactivity through a 
cross-beam scattering method [Slj, HI]. A particularly 
noteworthy contribution is the study by Shapiro and 
Brumcr of the coherent control of single-molecular pho- 
toassociation or bimolecular collisions through the inter- 
ference of reactive pathways [HI . 

Extending these considerations to the case of ultracold 
matter waves, we show that the coherent abstraction re- 
action can be realized efficiently and controlled in a STI- 
RAP photoassociation pulse sequence [l3| such that the 
intermediate states are dark states. An important char- 
acteristic of this process is that it is triggered by quantum 
noise, leading to large shot-to-shot quantum fluctuations 
that invalidate the use of the Gross-Pitaevskii equation 
(GPE) in the initial stages. That equation can be used 
only at later times when the product reactant channels 
become macroscopically occupied. This is somewhat sim- 
ilar to a situation familiar in a number of quantum and 
atom optics examples such as the laser [s^l, optical and 
matter- wave superradiance [ssj . and molecular matter- 
wave amplifier [361 . and is in contrast to the familiar 
single- molecular |lOl. [26l . [37j combination reaction. 

In realizing the collective reaction A + B2 ^ AB + B, 
the basic idea is to first create weakly bound trimers 
AB2 via an entrance-channel atom-dimcr FR, and then 
to dissociate them into a closed-channel bound dimer 
and atom via photodissociation. A key aspect of that 
scheme is that involving a trimer intermediate state 
allows one to exploit a coherent population trapping 
(CPT) state that prevents the trimer population from 
becoming significant throughout the conversion process. 
Such an atom-molecule state does not exist in other 
schemes that involve e.g. an intermediate two-species 
atomic state. Note also that this scheme, which is spe- 
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cific to quantum-degenerate matter waves, is different 
from a purely collision-induced reaction [3^ and from 
the non-degenerate single-pair dynamics of reactive scat- 
tering (33! . As such it represents a promising advance in 
the on-going development of superchemistry [3, HI] ■ 

The remainder of this paper is organized as follows: 
Section II discusses our model and analyzes the initial 
stages of the coherent bimolecular reaction. The result- 
ing quantum fluctuations determine the initial statisti- 
cal properties of the mean-field evolution that takes over 
once the various matter-wave fields are macroscopically 
populated. This is discussed in Sec. Ill, where we numer- 
ically compute a large number of trajectories from initial 
classical seeds satisfying these short-time statistics. We 
also review how an approximate CPT dark state required 
for the STIRAP pulse sequence can be achieved in the 
presence of mean-field shifts. Section IV briefly discusses 
the possible conversion of bosons to fcrmions. Several 
generalizations, including the role of population imbal- 
ance, are considered in Sec. V. Finally Section VI is a 
summary and conclusion. 



II. SHORT-TIME QUANTUM DYNAMICS 




FIG. 1: (Color online) Schematics of the coherent abstraction 
reaction A + B2 —> AB + B with degenerate matter waves. 
A and B denote the bosonic or fermionic atoms, and B2 and 
AB are molecular dimers. 



Our model system is sketched in Fig. 1. The inter- 
mediate hetcronuclear trimers AB2 are created via FR, 
and then photodissociated into bound molecules AB and 
atoms B. Denoting the atom-dimer coupling strength 
with detuning d by A'^, the Rabi frequency of the disso- 
ciating laser by fl'i and its detuning by A, the dynamics 
of the system is described at the simplest level by the 



model Hamiltonian (h = 1) 

+ A'[^4(r)^,(r)Vib,(r)+H.c.] -f (A + (5)^^,(r)46W 
- r!'[^L(r)V'l(r)V'™M+H.c.]}. (1) 

We consider first a purely bosonic system. The an- 
nihilation operators ipi{r), where the indices i,j — 
a,b,b2,ab,m stand for atoms {A and B), dimers {B2 
and AB) and trimers {AB2); satisfy the standard bosonic 
commutation relations 

[Mr), ^](r')]=%<5(r-r'). 

The terms proportional to Xi j = '^T^O'i.j / Mi^j describe 
interspecies s-wave collisions with scattering length a^j , 
Mi J = M,Mj/{M, + Mj) being the reduced mass 

In the mean-field approximation, ipi — > ^Jnil^i where n 
is the initial particle density, the Heisenberg equations 
of motion resulting from the Hamiltonian (1) are easily 
shown to reduce to the form 



■0a 


j 






j 






j 






j 


- inTplx/jjn + j(A + 5)lpab 


V'm 


= 2iExmj|V'iPV'm 


+ {iS - 7)'(/)„i + iXi>aipb2 



3 



-iflipbtpab- (2) 

Here Xi,j = "Xi,j' ^ = A'-^n, ^ = ^'\/n, and we have 
introduced the phenomenological decay rate 7 to account 
for the loss of intermediate trimers, based on the assump- 
tion that this decay dominates all other loss mechanisms 
such as rogue photodissociation to noncondensate modes 

|43|. As already mentioned, our goal is to minimize 
that decay by using a STIRAP pulse sequence, ideally 
permitting the full transfer of the entrance-channel state 
to the closed-channel state while keeping the intermedi- 
ate state unpopulated at all times. 

The initial condition ipb{0) ~ ipabiO) = is readily seen 
to result in 

Mt) = V'ah(t) = 0, 

for all times. This indicates that the mean-field GP 
equations break down completely in studying the onset 
of this type of abstraction reaction. A similar situation 
has been previously encountered in a broad range of sys- 
tems in quantum optics [3^ . but also in coupled degen- 
erate atomic and molecular systems such as the exam- 
ple in the matter-wave superradiancc of Bose-condenscd 
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atoms [3J, [33, [Sg]- As in those situations, our strat- 
egy here is to decompose the problem into an initial 
quantum-noise-dominated stage followed by a classical 
stage that arises once the product components have ac- 
quired a macroscopic population. The initial quantum 
evolution is treated in a linearized approach whose main 
purpose is to establish the statistical properties of the 
initial fields required for the classical stage [1^. 

A simple physical picture of the initial stages of the 
coherent abstraction reaction can be obtained when con- 
sidering the limiting case where S is the largest parameter 

in the system, iipm/S « 0. In the coUisionless limit this 
gives 

1pm « -{X'/S)lpa1pb^ + {n'/S)lpb1pab, 

which amounts to adiabatically eliminating the interme- 
diate trimer state. In this case the system is described 
by the effective Hamiltonian 



Hoff = -{Gclb^b^aCb2 + h.c.) + Co, 



(3) 



where 



Co = UJic\cacl_^Cb2 + ^2c\f,Cabc[cb-, (4) 



and the various constants are 



G ^ (X'n'/S) J dr<j,:,{r)rb{r)Mr)^b2{r). 
^1 = (A'V^) / drra{r)MrWbAr)<l>bAr), 
L02 = (n'yS) I drrab{r)4'abirWbir)Mr), 



The Hamiltonian (3), which is exactly solvable, has been 
considered previously in the study of a spin exchange 
scattering process that produces entangled bosonic pairs 
in a two-species, two-pseudospins Bose condensate |4l| . 

For short enough interaction times, Eq. (3) can be 
further simplified by taking into account the fact that 
the populations of the reaction products remain small 
compared to the total particle numbers A'o. In this 
regime, we can treat the fields "ipa and '>pb2 classically, 
Ca,b2 ~* \/Na,b2; ^ud thcn ncglcct the term in Eq. (3) de- 
scribing only the interactions between the modes and 
1)2- As discussed e.g. in Ref. [sl], this results in the lin- 
earization of the Heisenberg equations of motion for the 
remaining quantized matter- wave fields Cab and Cb, with 
a noise source fUt), 



'0ab,b{t)=flab{i)=M.ab{t)- 



(5) 



Such a result is familiar from several quantum and atom 
optics problems, including the optical parametric oscil- 
lator [3J| and molecular dissociation (pair production) 



in ultracold systems |42| . For simplicity we assume in 
the following that the correlations of the quantum noise 
operators appearing in Eq. ([5|) arc markovian, 



{fKt)fj{t')) - 0, {mn{t')) = gH,,5{t - 1') 



where Q = G^jNaWj^ and i,j = abovh here and in the 
following. It is these noise operators that trigger the non- 
mean- field "spontaneous" evolution of the system from 
initial vacuum fluctuations. 

The populations of the modes j — b, ah and their cor- 
relations are then 



Nj = (c]cj) = sinh' (gt) « NaNb2GH'' 



Cn.h — Cb 



^ {ANabANb) 



i + fiiYi\?{gt)>\, (6) 



where AiVj = Nj — (Nj). Eq. (5) can also be derived 
by solving Eq. (3) to second order in time t with the 
depletions iVa,b2= — A^ab,b- It is also straightforward 
to find that 



and that 
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It follows that the Mandel Q parameters [3J] is given by 



Qa 



b,b/ 



N 



b,b 



b.b 



N, 



= cosh^ (gt) > 1, (8) 



ab.b 



and exhibits super-Poisson statistics [3J| . It is interesting 
to observe that although the second factorial moments 
of the single modes ah and b are typical of chaotic fields, 
quantum entanglement within these two modes does ex- 
ist, i.e.. 



(2) 
'ab,b 



yab,b 



{NabNb) 
NabNb 



= 1 



cosh'^ {gt) 



(2) (2) 
9ab9b 



sinh^gt) 
sinh"'*(gt) +4sinh~^(gi) > 0, 



(9) 



violating the classical Cauchy-Schwarz inequality (CSI) 

Hi,!!!]. 

Similar equations of motion can be derived in case 
atoms A are bosonic and atoms B fcrmionic. The main 
difference in that case is in the commutation relations of 
the noise operators (— fib)- Heisenberg equa- 

tions of motion can be solved via a Bogoliubov trans- 
formation. One finds that the vacuum-noise-triggered 
populations of principal modes are then (4^ 

Nab.b = sin^ (gt) < 1, 

a direct consequence of the Fermi statistics [l^. Being 
similar to the bosonic case, we also find {Nab) ~ (-^b ) = 
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sin 

-sm^{gt)cos{2gt)), and 



sin2(5i), {NabNb) 



(2) (2) 
9ab = 9l 



2 

9ab.b 



(2) 



NabNb 



{Nab 



= 1 - 



0, 



cos^(gt) 

sin^{gt)' 



^.^fff = [1 - > 0. (10) 



The fermionic dimer-atom pairs correlations are of course 
also different from the bosonic case, specifically we have 
now Cab,b = 1 — Nab.b < 1, a signature of antibunching. 
Additionally, the Mandel Q parameter for the principal 
mode is 



Qab,b — 



{Nlb.b 



Nlbb 



= cos^ (gt) < 1, 



(11) 
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characteristic of subpoissonian statistics. 
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FIG. 2: (Color online) Standard derivation of the dimer and 
atomic populations from their average values for 5 = 3 and 
S — —3. Time is in units of A~^, and 7 = 1. The other pa- 
rameters are given in the text. The trimer population remains 
essentially zero at all times due to the CPT condition. Inset; 
fluctuating range of the populations Ni{t) ± A7Vi(t) for 5 = 3. 



III. LONG-TIME CLASSICAL EVOLUTION 

The long-time statistical properties of the AB and B 
populations, which arc significantly influenced by the ini- 
tial vacuum fluctuations, can be calculated by a positive- 
P representation technique [2^ and other methods [4]| . 
Rather than adopting such a full quantum treatment, we 
proceed in the following by solving the mean-field de- 
scription of Eqs. (2) with stochastic classical seeds whose 
statistics are consistent with the results of the linearized, 
short-time quantum analysis [1^. To be specific, using 
Eqs. (2) we compute 300 trajectories with randomly cho- 
sen initial classical seeds satisfying the short-time behav- 
ior of Eq. (5). 

Figure 2 shows the standard derivations AiV^ (i) around 
the average values of the particle populations, 

300 ~| 

with N,{t) = (1/300) J^n N^,n{t), for ,5 = 3 and 6 = -3. 
The inset shows the fluctuating range ztANi about the 
mean populations, NizL ANi{t) for 5 = 3 and for bosonic 
atoms. The small seeds resulting from the initial quan- 
tum fluctuations are significantly amplified, increasing 
more rapidly than their deviations, before reaching a sta- 
tionary value. For S = —3, however, no stable reaction is 
observed. 

An important feature of the coherent abstraction re- 
action is that it can be controlled and optimized by ex- 
ploiting the existence of a CPT dark state [1^, [IJ . This 
technique is well known in the case of linear systems, 
where it permits the transfer of population from an ini- 
tial to a final state via an intermediate state that remains 
unpopulated at all times. This is the basis for stimulated 



Raman adiabatic passage (STIRAP), which achieves this 

Pil via a so-called counter-intuitive sequence of pulses 

CPT and STIRAP rely explicitly on the validity of the 
adiabatic theorem, which applies only to linear systems. 
While there have been many previous developments in 
associative STIRAP (see e.g. Ref. [26]), it is not im- 
mediately obvious that STIRAP still works in the usual 
way with a well-defined nonlinear adiabatic condition for 
the nonlinear system at hand, even in the coUisionless 
limit. This question was recently investigated by Pu and 
coworkers in Rcf. [16] for a model of coherent atom-dimer 
conversion, and an approximate adiabatic condition was 
obtained by linearizing the nonlinear system around the 
intended adiabatic evolution. In that context, the ana- 
lytical form of the nonlinear adiabatic condition derived 
in the coUisionless limit turns out to be useful in deter- 
mining the required laser parameters [l^, [13, d. El] . 

We show in the following that an approximate atom- 
moleculc CPT state can also be achieved in the present 
situation. Specifically, Eqs. (2) admit a steady-state CPT 
solution with a trimer state that remains unpopulated at 
all times under the generalized "two-photon" resonance 
condition 

A = -S + {2Xaa + GXab + SXfcbs + '^Xbb)Nb2,s 

+ {2Xab + Xbh2)Nab,s, (12) 

To this end, we apply the steady-state ansatz 

^a,s = ]^a,sle''"e^''°*, i^b,s = \i^b,sW''e'''''. 

i>b,.s = \i^b,Ae^'''e^'^'\ 

i^ab.s = ]V^ab,.]e^('"+'^)e^('^"+^'')*, 

^™,. = lV™,.le^('"+^'^^)e^('^"+2^'')*, (13) 
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TABLE I: Steady-state CPT particle numbers, with Sle// 



where /i^ and /it are the atomic chemical potentials. In- 
serting these trial functions into Eqs. (2) and taking 
|V'm,s| = 0, one finds the steady-state relation 

which, together with the condition of conserved particle 
numbers Na,s + 2(iVb2,s + Nab,s) + A^fc,.s = 1, gives 



where 



Na,s + Nab.s 



The CPT solution is therefore 

27^ 



= 0, (14) 



(15) 



N, 



ab.b 



(1 + 7^) [i + 27e + vTi~^2^F+8mF7A2] 



Table 1 displays the CPT particle numbers for several 
values of TZ. Note that from the conservation of particle 
numbers we have Nb s = Nab.s and 



N, 



1 



2{i + n) 



In addition, we note that an initial populations imbal- 
ance can also significantly affect the dynamics of atom- 
molecule conversion for our present bimolecular reac- 
tions. The condition d Nab.s /dTZ =^ yields a maximum 
dimcr number 



Nab.s 



1/3, 



corresponding to a complete abstraction reaction for 
TZ = 1/2 or the so-called "balanced case". We consider 
this case first and then turn to study the effect of initial 
populations imbalance (see Section V). 

We have numerically solved Eqs. (2) and the typical 
results are showed in Fig. 3. In this specific example 
atom A is ®^Rb, atom B is -^^K, A = 4.718 x lO^'s"! and 

n{t) ^ f}osech(t/r), 

with r^o/A = 20 and At=20 [sO]. The collision param- 
eters, in units of X/n, are Xaa = 0.5303, Xbb = 0.3214, 
Xab = 0.8731, aU others being equal to 0.0938 As 
mentioned earlier, we have neglected rogue photodissoci- 
ation to noncondcnsate modes, which is proved through 



our direct calculations to be a safe approximation for 
the present parameters of our model |5ll . [s^ . [ssj . On 
the other hand, we note that the scattering lengths of 
the various particles collisions, especially those involving 
molecular dimers or even some trimers, depend on the 
details of the interatomic potential, are yet not known. 
However, in our numerical calculations it is straightfor- 
ward to use a large set of plausible collision parameters 
for the Rb-K, Rb-Na or other alkali atomic samples. We 
actually have done this and found essentially the similar 
result as Fig. 3: stable bimolecular conversion is always 
possible for appropriate values of the external field detun- 
ing (5, which is independent of the precise collision values. 
This result finds its origin in the underlying mechanism of 
"generalized STIRAP,' first proposed by Ling et al. [lB| 
in the context of atom-dimer conversion, and according 
to which collisions need not limit the conversion rate as 
long as one chooses an adiabatic passage route that com- 
pensates for the collisional mean-field phase shifts (see 
Eq. (12)). 

Figure 3 shows the creation of AB and B for i5 = ±3 
and (5 = ±1. The results are essentially the same as those 
of Ref. [s^] for i5 = ±3. For this bosonic system, sta- 
ble bimolecular conversion is always possible for negative 
detunings, but the system can be unstable for positive 
detunings (see Fig. 3). The increasing departure of the 
product populations from the ideal CPT line is due to the 
fact that only an approximate adiabatic condition exists 
for the CPT state. 

We also analytically derived the adiabaticity param- 
eter introduced by Pu et al. in the coUisionless limit 

Em, 



7nl(t) 



1^1 



1 



1 -f 77 4A 



< 1, 



where ry = A/f2. This expression differs from that for the 
corresponding linear system [l6| in that in the latter case 
77 is replaced by 77^ in the denominator. Hence adiabatic- 
ity becomes increasingly difficult to maintain in the final 
stages of the STIRAP process, similarly to the case of 
atom-dimer conversion jl6| . 

We remark that this scheme relics crucially on the ca- 
pability to avoid rapid collisional quenching or the forma- 
tion of an unstable atom-dimer sample. When energeti- 
cally allowed, collision-induced reactions always occur at 
some rate, and we need to guarantee that the time scale 
over which quantum fluctuations dominate the dynamics 
is short enough, so that the dynamics of the system is 
not collision-dominated. 

In order to estimate the upper limit on collisions, we 
can use the condition \Qt\ < 1 or 



IGy/NM^tl < 1, 

which determines the validity of the short-time approx- 
imation for the early quantum stage, to estimate the 
time over which the fluctuations take place. This up- 
per time limit is of the order of 10~^s for |(5 = 3| and 
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FIG. 3: (Color online) Populations of the bosonic dimers and 
atoms for (a) & — ±3 or (b) S — ±1. Time is in units of 
A~\ and 7 = 1. The line labelled "CPT" shows the ideal 
population of products (dimers AB and atoms B). 



flo = 20 A. According to Cvitas et al. [5^, typical low- 
temperature inelastic collision cross-sections are of the 
order of 10~^^m'^/s, corresponding to reaction times of 
the order of 10~''s for a typical condensate density of 
lO^'^cm"'^. From that estimate, it appears that the fluc- 
tuations do indeed dominate for short enough times. We 
also note that the coUisional reaction time of lO^'^s corre- 
sponds to an almost complete noise-amplified conversion 
in Fig. 3. In that case the fluctuation-induced dynam- 
ics completely dominate the short-time behavior of the 
system. This feature, which is characteristic of a wide 
variety of collective abstraction reactions, may provide a 
useful means to produce reaction products that are dif- 
ficult to obtain or have only poor yield when resulting 
from a purely coUisional method. 



IV. BOSONS TO FERMIONS CONVERSION 

When considering a mixture of bosonic and fcrmionic 
atoms, the abstraction reaction results in the conversion 
of bosonic to fcrmionic molecules, 

where b {B) or f {T) denotes bosonic or fcrmionic atoms 
(dimers). Ignoring s-wave collisions between fcrmionic 
atoms of the same species and retaining only their domi- 
nant kinetic energy, and assuming further that collisions 
arc the dominating term for the bosons (soj , the Hartree 
energy density of the system is 

+ X' [^di^ai^b^ + h.c] ~n' [i^lbiJl^d + h.c] 



(16) 



i—{a.b2 -d) 



f={ab,b) 



Here the indices i, j have the same meaning as in Eq. (1) 
and A'j = (67r2)2/3/2Af/, where Mf (/ = ah,b) is the 
mass of the fcrmionic components. 

Due to the fcrmionic components of the reaction part- 
ners, the ordinary mean-field approach [s^, HI] adopted 
in studying systems with large numbers of condensed 
bosonic particles is inadequate here. We follow instead 
the approach of Ref. [5^, starting from the mean-field 
Lagrangian density of the system 

- = ^E(*.-f-*f)- m 

i 

and exploiting the Euler-Lagrange equations 

1=0. (18) 



9£ / dC 

9V* ''va(a^V'r) 



to derive the mean-field dynamical equations 

3 

tjjb =2i"^xbj\i^j\^i^b + iAbliJbl'^/^'ipb - li^iPlbi^d, 

i 

i^ab =2i ^ XafcjIV'jI^V'oh + iAab\tpabf^^1pab 
j^ab 

il'd=2i^Xd,j\tpj\'^ip,n + {iS - 'y)i',n + iXtl'atpb2- (19) 
j 

These equations are similar to Eqs. (2) with the substi- 
tution 



(20) 
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a consequence of the fact that we consider only the dom- 
inating kinetic energy term ignore s-wave coUisions be- 
tween identical fermionic particles [sHi ■ 

In the CPT regime, the steady-state number of 
fermionic species AB and B is therefore in the same 
form as in the purely bosonic case, see Table I. How- 
ever, the generalized "two-photon" resonance condition 
is now written as 

A = -S + 2(Xab + Xaa + Xab2)^b2,s + 4XabiVab,s 

+ (Ab - ^ab)A^f£. (21) 



We note that in contrast to the purely bosonic case, 
the creation of fermion-fermion pairs is due to a statistics- 
independent cooperating many-body effect that has been 
previously recognized in the case of, e.g., matter- wave 
four- wave mixing [s^. We also note that in recent work 
Li et al. (60| used a similar atom-molecule dark-state 
technique to realize a so-called laser-catalyzed bimolec- 
ular reaction (or the conversion of fermionic to bosonic 
molecules) : 



Li +^ Li^Li Li2 -1-^ Li, 



0.4 



(a) Boson A and Fermion B 



CTP (AB, B) 




AB^B (6=-3) 
AB7 6(5=3)" 



A, B^ (6=3) 



50 100 150 200 250 300 
time 



(b) Boson A and Fermion B 



CPT (AB,B) 




AB, B (5=-1) 
AB~b"(5=1)~ 



A, B^ (5=1) 



50 100 150 200 250 300 
time 

FIG. 4: (Color online) The generation of fermionic dimers AB 
and atoms B for (a) 5 — ±3 or (b) 5 = ±1 initially with no 
particle populations imbalance. The CPT values of the dimer 
AB and atomic B are also plotted. 

Our numerical simulations of this case are summa- 
rized in Fig. 4 for the collision parameters Xaa = 0.5303, 

Xbb = 0, Xab = Xb,ab = Xa,ab = -0.09, Xa,ab = "0.2637, 

Ab = 0.008, Aab = 0.004, and all other collision param- 
eters equal to zero. These results show that the stable 
formation of dimers AB is possible for both positive and 
negative detunings ((a)(5 = ±3 or {h)S = ±1), initially 
with no particle populations imbalance. 



with an ultrahigh conversion rate of 99.97% [6C 



ROLE OF POPULATION IMBALANCE 



Population imbalance often plays an important role in 
the physics of ultracold matter waves. For example, us- 
ing a two-spin-state mixture of ultracold fermionic atoms, 
population imbalance can induce a supcrfluid to normal 
phase transition [6l| . In the following we demonstrate 
that an initial population imbalance can also significantly 
affect the dynamics of this coherent collective abstrac- 
tion reaction. To this end we plot the conversion rate as 
a function of TZ, see Eq. (14). Section III showed that 
in the collisionless limit Nab,s reaches its maximum for 
TZ = 0.5. This value is modified slightly when taking into 
account of the particle collisions and the decay rate 7 in 
our numerical simulations. 

In Fig. 5, we see that the conversion rate \tpabit = oo)p 
now has the maximum at the value somewhat larger than 
R — 0.5 for the both bosonic and fermionic cases. In 
addition, we observe that the initial population imbal- 
ance has different effects for bosons and for fermions: for 
bosonic atoms A and B, the final conversion rates can 
be changed quite sharply with IZ and rapidly approach 
zero for TZ < 0.4 and TZ > 1.4; in contrast, for the case of 
fermionic atoms B, the conversion rate changes its shape 
more slowly with TZ (even for TZ > 1.4, the occupation of 
the product species is still in excess of 7%). 

Finally, we remark that we can follow a similar ap- 
proach to study the collective multi-molecular reactions 
2AB A2 + B2 and 2A2 A^i + A. In the first case 
we find that the steady-state CPT values of the prod- 
uct dimer A2 or B2 are the same as that of AB in the 
A + B2-^ AB + B reaction for 7^=1/2, see Table 1. For 
the reaction 2A2 — * A3 -I- A, we obtain the equations of 
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tion as 
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(b) Boson A and Fermion B 
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FIG. 5: (Color online) Final dimer population \'^ab{t = c>o)|'^ 
as the function of R for (a) the Bose case and (b) the Fermi 
case for several values of 5. 



motion 

3 
3 

V'aa ='^i'^Xa3,3\^j?^o.3 + i(A + 5)?/'a3 ^ i^'4^*a'^t 
3 



(22) 



from which it is straightforward to show that the "two- 
photon" resonance condition is 

A = -,5 +(3Xa.,a " 2Xa,,a)|V'al' + (3Xa.,a. " '^Xt.a.Waf 
+ {-iXa,,a,-2xt.a,)\raf. (23) 

Introducing a steady-state ansatz similar to Eqs. (12), we 
can then calculate the CPT value of the trimer popula- 



1 

2(2 + r2/A)' 



(24) 



We have carried out numerical simulations for this case 
and found results similar to those sumarized in Fig. 3. 



VI. CONCLUSION 

In conclusion, we have shown that in quantum- 
degenerate ultracold atomic and molecular systems the 
abstraction reaction A + B2 ^ AB + B can be funda- 
mentally different from the familiar atom-molecule com- 
bination. The quantum noise that triggers this reaction 
leads to the creation of strongly correlated dimer-atom 
pairs [4l| . In addition, a generalized atom- molecule dark 
state existing in this system can significantly enhance the 
creation of boson-boson or fermi-fermi dimer-atom pairs. 
Our study opens up a new direction in association reso- 
nances and thus a fascinating new regime of superchem- 
istry. In particular it can be generalized to the analysis of 
more multi-molecular matter-wave reactions, such as the 
reaction 2AB A2 + B2 or the creation of a quantum- 
entangled atom-molecule laser A + A + A ^ A + A2 [12, 
39, 62], controlled even by purely optical means. 

Future work will investigate the unique "superchem- 
istry" effects of ultra-selectivity or confinement-induced 
stability on the bimolecular reactions, the collective ab- 
straction reaction in a double- well potential or an optical 
lattice, and the possible applications of these kinds of re- 
actions in realizing laser-catalyzed atomic spin mixing in 
a spinor-1 Bose condensate |27l |63|. In addition, a com- 
plete analysis of coUisional effects |38l. l64j will be consid- 
ered. While experiments along the line of this analysis 
promise to be challenging, recent progress in quantum 
degenerate chemistry jl5» .25,, 211 a^nd in the manipula- 
tion of atom-molecule systems [l^, [13, HO, HI, Hi, 
indicates that achieving this goal should become possible 
in the not too distant future. 
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